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The ensemble of L x L power-law random banded matrices, where the random hopping Hij 
decays as a power-law {b/\i — jl)", is known to present an Anderson localization transition at a = 1, 
where one-particle eigenfunctions are multifractal. Here we study numerically, at this critical point, 
the statistical properties of the transmission T2 for two distinguishable particles, two bosons or two 
fermions. We find that the statistics of T2 is multifractal, i.e. the probability to have T2{L) ~ l/L'' 
On ' behaves as where the multifractal spectrum $2(k) for fermions is different from the common 

' multifractal spectrum concerning distinguishable particles and bosons. However in the three cases, 

, the typical transmission Tn^^{L) is governed by the same exponent ^2^^, which is much smaller than 

the naive expectation 2^-^"^, where k'^^ is the typical exponent of the one-particle transmission 

Oh! 
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I. INTRODUCTION 



Whereas Anderson localization phenomena [l| are rather well understood for a single particle (see the reviews 
rt^' ^^^^ of interactingparticles in a random potential has remained much more challenging (see the review 
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3\ and more recent works [gI. Isl. loL [loL [llj| ) . Since the case of a finite density of particles can be studied numerically 
only for small system sizes, it is natural to consider first the simpler case of only two interacting particles (T.I. P.) 
in a random potential. In dimension d = 1, where the one-particle model is always in the localized phase with some 
localization length Ai, it has been found that the T.I. P. is also always localized, but with a localization length A2 
that may become much larger than Ai [H, [H, [3, EH [3 [I3j [H) IS 113) HH- In dimension d = 2, where the one- 
particle model is again always in the localized phase, the possibility of a delocalization transition has been studied 
for short-range interaction [2§| and for Coulomb interaction [l^, [23| ■ 

In the present paper, we are interested in the two-particle transport properties at the Anderson localization tran- 
sition of the one-particle problem, where the one-particle eigenstates arc multifractal 0, We are not aware of 



O ' previous studies on this question (see however 25| concerning the quasi-periodic Aubry- Andre transition). Since for 

^ O ^ I the tight-binding model in dimension d = 3 where there exists an Anderson transition, the two-particle model cannot 

I be studied numerically for large enough system sizes and large enough statistics on the disordered samples to obtain 

^— I . accurate results, we have chosen to focus here on the Power-law random banded matrices (PRBM) model, and to 

J> ' study numerically the statistical properties of the two-particle transmission T2. 

, The paper is organized as follows. In section |TT1 we recall the Power-law random banded matrices (PRBM) model 
and introduce the obscrvablcs that characterize transport properties for the two-particle model. In section IIIIl we 

••^ , describe our numerical results concerning the statistical properties of the transmission T2 for two distinguishable 

. • particles, two bosons and two fermions. Our conclusions are summarized in section ITVl In Appendix [Xl we describe 

On ' our numerical results concerning the multifractal properties of the one-particle model as a function of the energy E, 



which turn out to be useful to understand the statistics of T2 discussed in the text. 



II. MODEL AND OBSERVABLES 



5_( ' A. Reminder on the Power-law random banded matrices (PRBM) model 

a ^ 

Beside the usual short-range Anderson tight-binding model in finite dimension d, other models displaying Anderson 
localization have been studied, in particular the Power-law Random Banded Matrix (PRBM) model, which can be 
viewed as a one-dimensional model with long-ranged random hopping decaying as a power-law (5/r)° of the distance 
r with exponent a and parameter b. The Anderson transition at a = 1 between localized (a > 1) and extended 
(a < 1) states has been characterized in (26| via a mapping onto a non- linear sigma- model. The prop erties of the 
critical points at a = 1 have been then much studied, in particular the statistics of eigenvalues [27l. l28l. [29j . and the 
multifractality of eigenfunctions [13, HU, i32, 33, 34, 35], including boundary multifractality [36| . 

More precisely, we consider here the model shown on FiglU with L sites i = 1, 2, ..L in a ring geometry with periodic 
boundary conditions. The appropriate distance r^.j between the sites i and j is defined as [30| 

(L) L ■ ( 7r(i - j)\ 
<,--sin(— ^) (1) 
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FIG. 1; (Color on line) The ensemble of power-law random banded matrices of size L x L can be represented as a ring of L 
sites, where the matrix element Hij between the sites i and j is a Gaussian variable of zero-mean Hij = and of variance 
given by Eq. [5] in terms of the distance rij of Eq. [T] 



The ensemble of power-law random banded matrices of size L x L is then defined as follows : the matrix elements 
i/j^^^ are independent Gaussian variables of zero-mean = and of variance 

The most important properties of this model are the following. The value of the exponent a determines the localization 
properties [26] : for a > 1 states are localized with i nteg rable power-law tails, whereas for a < 1 states are delocalized. 
At criticality a—1, states become multifractal (sol l3ll [s^. [ssj and exponents depend continuously of the parameter 
6, which plays a role analog to the dimension d in short-range Anderson transitions [s^l : the limit 1 corresponds 
to weak multifractality ( analogous to the case d = 2 + e) and can be studied via the mapping onto a non-linear 
sigma-model (26|], whereas the case 6^1 corresponds to strong multifractality ( analogous to the case of high 
dimension d) and can be studied via Levitov renormalization [sol . Is^ . Other values of b have been studied numerically 
[sol , [sil . [3^ . [3^ . The statistical properties of the Landauer transmission for a single particle between the opposite 
points L/2 and L has been studied in detail in our previous work i^Sj (results concerning other scattering geometries 
can be found in [39j). 



B. Transmission of two distinguishable particles, two bosons, or two fermions 

In this paper, we consider the two-particle model defined by the Hamiltonian 

i7(2) = ii-(i) 1 + 1 (3) 

As stressed in [H |T3, [H, ^ [2l| for the one-dimensional T.I. P. model, the important observable to characterize the 
two-particle transport properties is the Green function 

between doubly occupied sites along the diagonal ri = r2. In our present notations concerning the P.R.B.M. model 
(see Fig. [Ij, we will thus focus on the transmission 

T,^\<^,^\Ge.^o\L,L>\^ (5) 

at zero energy £'2 = (center of the band). It is important to stress that even if there is no explicit interaction in the 
Hamiltonian of Eq. [31 the two-particle Green function cannot be factorized into one-particle properties [13, HO] ■ We 
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will indeed find below non-trivial properties for T2. As a comparison, one may also consider the transmission of one 
of the two particles with the other held fixed (see Eq. 7 of |19j ) 

r2,(i/)^|<^,i|G|i,i>|' (6) 



III. NUMERICAL RESULTS ON THE STATISTICAL PROPERTIES OF Ta 



A. Numerical procedure 

We have used an exact diagonalization method the one-particle P.R.B.M. model H^^^ for the localization transition 
critical value a = 1 and for the parameter b = 0.1 (Eq. [J). For each disordered sample, we note e„ the L eigenenergies 
(n = 1, 2, ..L) and the corresponding normalized eigenstates {i — 1, 2, ..L) 

H'^^^ni^ = e„0„(^) (7) 

To compute the two-particle Green function, we now have to know the symmetry properties with respect to the 
exchange of the two particles. 



1. Two distinguishable particles (No symmetry conditions) 

For two distinguishable particles, an orthonormal basis of eigenstates of H^"^^ is given by the following states 
labelled by two integers \ <n < L and 1 < m < L 

■fpn,m{i,j) = (l)n{i)(j)m{,j) (8) 

of energy 

En^m — ~t- C^n (9) 

The two-particle Green function at zero energy £'2=0 then reads 

G'£;,=o(«,j;« ,j ) =-2^2^ p (10) 



n— 1 rn—1 
L L 



n.7n 



J2 ^ Km*mU)Mi')^m{f) ^^^^ 



2. Two Bosons (Symmetry condition) 

An orthonormal basis of eigenstates is given by the following L{L + l)/2 symmetric states labelled by two integers 

1 < n < m < L 

<n(*,i) ^M^)Mj) (12) 

iB r ■\ _ 4'n{i)4>ni{j) + 4> 0') 

of energy given by Eq. [HI 

The two-boson Green function at zero energy £'2=0 then reads 

G£,=o(*,j;«j) = -2^ 2^ -— (14) 

1 \ 
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3. Two Fermions (Antisymmetry condition) 

An orthonormal basis of eigenstates is given by the following L{L — l)/2 antisymmetric states labelled by two 
integers 1 < n < m < L 

Wn,m.V'^3) — U^j 

of energy given by Eq. [5] 

The two-fermion Green function at zero energy E2 = then reads 

I- ■ -t ■t\ V'n,m(*) j)V'n,m(* ) 

G£,=o(*-j;m ) = -2^ 2^ -— (16) 

T , T tin ~r Cm 

n— 1 m— n+1 

For fermions where double occupancy is forbidden, we have modified the definitions of Eqs [S] and [S] for the trans- 
missions into 

= |<L/2,L/2-l|G|L,L-l>|2 (17) 

and 

T2^(i^)^|<L/2-l,L|G|i:-l,L>p (18) 

The results given below correspond to sizes 50 < i < 2000, with corresponding statistics of 5.10^ > ns{L) > 1150 
independent samples. To improve the statistics, we have considered, for each disordered sample, the transmission T2 
between the L/2 pairs of opposite points. All results concern the zero-energy {E2 = 0) transmission T2 at the critical 
point a — 1 and the value & = 0.1 (see Eq. [2]). We first focus on the scaling of the typical transmission before we turn 
to the multifractal spectrum. 

B. Typical transmission T2^^{L) as a function of L 

We find that the typical two-particle transmission 



T2*^^'(L) = e'"^^(^) (19) 

decays as the power-law 

Tr{L) ^ (20) 

where 

~ 1.86 (21) 

is the same for two distinguishable particles, two bosons or two fermions as shown on Fig. O 

As a comparison, we also show on Fig. [5] the typical transmission T^^ij) of El- [HI representing the transmission of 
one of the two particles with the other held fixed : for distinguishable particles, bosons or fermions, it is governed by 
the same exponent 



K. 



Jyp 

2,(1/) 



1.3 (22) 



that coincides, within our error bars, with the exponent measured in [38l | for the one-particle model. 

C. Multifractal statistics of T2 

We find that the statistics of T2 is multifractal, i.e. that the probability to have T2{L) ~ l/i" behaves as 

Proh (T2{L) ^ L-^) dT cx L'^^^'^Uk (23) 



5 




FIG. 2: (Color on line) Scaling of the typical two-particle transmission T2^^ at criticality a — 1 for b = 0.1 and at zero-energy 
E ^ : \nT*y^{L) = In7^(I) as a function of InL for fermions (□), bosons (O) S'lid distinguishable particles (x) yield the 
same exponent k^J^'' ~ 1.86 (Eq. I2ip . As a comparison, the scaling of the typical transmission T^^^^^ of Eq. [S] (representing 

the transmission of one of the two particles with the other held fixed) corresponds to the exponent Kj^i/) — 1-3 ( Eq. [22}. 




FIG. 3: (Color on line) Multifractal spectra ^2(1^) describing the statistics of the transmission T2 at criticality a — 1 for fe = 0.1 
and at zero-energy E — : for two fermions (□), two bosons (Q) and two distinguishable particles (x). As a comparison, the 
multifractal spectrum $i(k) describing the statistics of the corresponding one-particle transmission is shown (o). 
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We show on Fig. [3] the multifractal spectra ^2{k) corresponding to two distinguishable particles, two bosons and 
two fermions. We find that the spectra for distinguishable particles and bosons coincide, whereas the spectrum for 
fermions is clearly distinct, except around the maximum ^2(1^*2^) = associated to the same typical value k^2^ of 
Eq. [21] This can be explained as follows : the transmission T2 for distinguishable particles and bosons both involve 
coinciding points (Eq. [5]), whereas the transmission T2 for fermions involves neighboring points (Eq. I17p . Besides 
their common typical scaling, one thus expects differences in their statistics. 

As a comparison, we also show on Fig. [3] the multifractal spectrum ^i{k) describing the statistics of the corre- 
sponding one-particle transmission Ti. A natural question is of course whether the multifractal spectrum $2(k) can 
be related to $i(k) or to the singularity spectrum of one-particle eigenfunctions. 



D. Discussion : relation with the statistics of one- particle eigenfunctions 



We first recall the case of the one-particle transmission, before we turn to the analysis of T2 



1. Analysis of the one-particle transmission in terms of one-particle eigenfunctions 



In terms of the energies e„ and eigenfunctions of the one-particle model (Eq. [7]) , the one-particle zero-energy 
Green function reads 



(24) 



In the limit of large size L where the levels become dense, the zero-energy Green function of Eq. [TT] become 

gE,=o{i,i)^-L dep{e) 



which is dominated by the neighborhood of e = 
so that the one-particle transmission scales as 



(25) 



(26) 



(27) 



When the distance \i — i'\ is of the order of the system size L, the weights |0e=o(*)P and |0e=o(*')P can be considered 
as independent. Then the multifractal spectrum $i(k) describing the distribution of the one-point transmission 



can be written as (here with d = 1) 



Prob (Ti - L-'") dT cx L'^^^'^Uk 



^i{K>Q) = 2\f{a = d+-)~d 



(28) 



(29) 



in terms of the singularity spectrum /(a) of zero-energy eigenfunctions (see more details in [38l. |40|). 



2. Analysis of the two-particle transmission in terms of one-particle eigenfunctions 



We now try to analyse the two-point transmission T2 for two distinguishable particles along the same lines. In the 
limit of large size L where the levels become dense, the zero-energy Green function of Eq. 1111 become 

Ge2=o(«, «; » , » ) dep[e) de p{e ) ; — (30) 



which is dominated by the region e + e' — 

Ge,=o{i; i) ~ -L'" I dep2(e)</):(z)Ce(«)0e(«')0-e(«') (31) 
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(where we have used the symmetry p(e) = p(— e) around the center of the band e = for the one-particle density of 
states). The two-particle transmission is then expected to scale as 

r2(z,z') = |Gs,.o(^;^')l' - L^" j dep\e) J de' p\e') [r,{^)r_,{^)^*A^W_A^)] x [</>e(*')0-e(*')0e' (*')0-e' (*')] (32) 

So we do not expect any simple expression for the multifractal spectrum $2 (n) : firstly, T2 contains eigenfunctions of 
any energy e, and the singularity spectrum /(a) of one-particle eigenfunctions depends continuously on the energy e 
(see more details in Appendix |A| ; secondly, T2 involves complicated correlations of eigenfunctions of various energies 
(studies of two-eigenfunctions correlations can be found in i42i| ) . 

It is however natural to consider the simplest approximation : if the integrals in Eq. [32] were dominated by 
e = — e' , one would obtain a direct relation with the one-particle transmission of Eq. [27] 

j.approx(a)(^^^,^ ^ L^'^ 1 0e=O (0 H </>e=0 (*') l' ~ iTl{^,^')? (33) 

In particular, the typical exponent k*2^ would read 

n'yP = 2K{yP (34) 

Our numerical results described above fEqs ?IT\ and [^^ show that k^^'' is in fact much smaller than (2k**'^). Our 
conclusion is thus that this simple approximation is very bad, and that correlations between one-particle eigenfunctions 
at various energies play a major role in the two-particle transmission T2. 



IV. CONCLUSION 



In this paper, we have studied numerically the statistical properties of the two-particle transmission T2{L) at the 
critical point of the PRBM model where one-particle eigenfunctions are known to be multifractal. Our conclusion 
is that 22 (L) is multifractal i.e. the probability to have 22 (L) ~ l/L" behaves as L*2(k)^ where the multifractal 
spectrum $2('*) for fermions is different from the common multifractal spectrum concerning distinguishable particles 
and bosons, because the double occupancy of a single site and the occupancy of two neighboring sites have different 
statistics at criticality. However in the three cases, the typical transmission T2^^(L) is governed by the same exponent 
K*2'^, which is much smaller than the naive expectation 2k**'^, where k^^^ is the typical exponent of the one-particle 
transmission Ti{L). This suggests that T2{L) probes non-trivial correlations of one-particle eigenfunctions of various 
energies. 



APPENDIX A: MULTIFRACTAL STATISTICS OF THE ONE-PARTICLE TRANSMISSION AS A 

FUNCTION OF THE ENERGY 

As first discussed in [4Qi] for the special case of the two dimensional quantum Hall transition, the critical probability 
distribution of the one-particle transmission Ti at an Anderson transition critical point takes the form of Eq. [28l 
where the multifractal spectrum ^i{k) can be related to the singularity spectrum f{a) of critical eigenstates via Eq. 

m 

For the P.R.B.M. model, numerical results on ^i{k) can be found in [38| at the critical point a = 1 and at zero 
energy e = for various values of the parameter b (Eq. [2]). 

Here we show on Fig. [J] how the multifractal spectrum ^i{k) at criticality a = 1 for the value & = 0.1 changes as a 
function of the energy e. In particular, the corresponding typical values read 



tif^ie = 0) 


~ 1.33 




P(e = 0.5) 


~ 1.38 




Kf^ie = 1) 


~ 1.55 




Kl^^ie = 2) 


~ 1.89 


(Al) 



Via Eq. [29l this shows that the singularity spectrum f{a) of critical eigenfunctions changes with the energy e. (The 
dependence on e of f{a) has been studied in [43|] for quantum Hall wavefunctions as a function of the Landau level). 
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FIG. 4: (Color on line) The multifractal spectrum "1?i(k) describing the statistics of the one-particle transmission (Eq. I28|l at 
criticality a = 1 and b = 0.1 for four values of the energy : e = (0)i e = 0.5 (□), e = 1 (o) and e = 2 (A). 

Since the zero-energy two-particle transmission T2 of Eq. l32l contains one-particle eigenfunctions of various energies, 
that are characterized by different multifractal singularity spectra, we do not expect any simple expression for the 
multifractal spectrum $2(k) of T2. 
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